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ABSTRACT 


A  mathematical  model  for  predicting  the  motions  of  SWATH  ships  in 
waves  in  the  frequency  domain  is  presented.  This  mathematical  model  requires 
definition  of  only  geometric  and  mass  properties.  While  the  mathematical 
model  for  the  vertical  plane  presented  here  differs  little  from  the  one  which  has 
been  used  for  about  a  decade,  the  transverse  plane  model  differs  from  the 
previous  model  in  various  ways.  Predictions  are  compared  to  measured  model 
experimental  responses  to  regular  waves.  These  results  show  that  the  responses 
to  regular  waves  are  well-predicted  for  a  range  of  hull  forms  and  conditions. 

ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  a  range  of  sponsors  over  many  years.  Most  recently,  the  effort  was  funded  by  the 
Naval  Sea  System  Command  (NAVSEA)  under  the  direction  of  the  Carderock  Division,  Naval  Surface  Warfare 
Center  (NSWCCD)  SWATH  Program  Office,  Code  1237,  under  Program  Element  64567N,  Task  Number 
AA10357,  and  Work  Unit  Number  1-1237-326. 

INTRODUCTION 

A  mathematical  model  for  the  motions  of  a  SWATH  ship  in  waves  has  been  developed  over  a  period  of 
years.  As  described  in  reference  1,  the  mathematical  model  used  for  SWATH  responses  to  waves  follows  the  strip 
theory  of  Salvesen,  Tuck,  and  Faltinsen^.  Lee^  applied  this  theory  to  twin  hull  configurations  and  utilized  an 
expression  presented  by  Thwaites"^  to  develop  a  model  for  the  cross  flow  drag  and  body  lift  contributions  to  the 
forces.  Hong^  introduced  surge  into  the  model.  Subsequently,  McCreight  and  Stahl^  developed  semi-empirical 
expressions  for  the  cross  flow  drag  and  lift  contributions  for  the  vertical  plane  responses  and  added  the  effect  of 
downwash  on  the  lift  of  the  aft  stabilizers.  Accurate  modeling  of  the  vertical  plane  motions  was  the  focus  of  this  last 
effort;  no  changes  to  Lee’s^  mathematical  model  for  the  transverse  plane  motions  were  made.  Recently,  the 
mathematical  model  for  transverse  plane  motions  and  the  appendages  has  been  modified.  Damping  terms  now 
couple  heave,  pitch,  and  roll. 

In  this  report,  the  six  degrees  of  freedom  mathematical  model  for  the  motions  of  SWATH  ships  in  waves 
which  is  implemented  in  SWMP94  (SWATH  Motions  Program,  1994  Version)  is  described.  With  this  model,  only 
geometric  and  mass  properties  are  required  in  order  to  predict  the  six  degrees-of-freedom  motions  of  SWATH  ships. 
Comparisons  between  model  scale  experimental  results  and  predictions  are  presented  for  vertical  and  transverse 
plane  motions  for  various  hull  forms. 

BACKGROUND 

Whereas  Lee^  and  Hong^  use  the  Frank  Close  Fit  Technique^  to  evaluate  the  velocity  potentials,  the 
computer  software  described  and  utilized  by  McCreight^  and  McCreight  and  Stahl^  includes  approximations. 
Velocity  potentials  are  not  evaluated.  Instead,  added  mass  and  damping  coefficients  are  approximated  using 


1 


expressions  developed  and  reported  by  John  Dalzell*.  Then,  with  various  assumptions  and  approximations,  the 
wave  exciting  forces  and  moments  are  approximated  as  a  function  of  added  mass  and  damping.  A  derivation  for 
approximate  heave  and  pitch  exciting  forces  is  given  in  reference  6,  but  the  derivation  for  the  transverse  plane  forces 
which  has  been  in  use  is  not  available. 

The  approximate  approach  for  the  added  mass  and  damping  coefficients  of  SWATH  configurations  was 
developed  by  Dalzell  in  the  late  1970's.  This  approach  was  advantageous  because  it  significantly  reduced  computer 
time  (and  cost)  compared  to  Frank  Close  Fit  Technique  calculations.  Consequently,  it  facilitated  evaluation  of 
numerous  hull  forms  in  design  studies.  The  approximations  assumed  that  the  cross  sections  had  wall-sided  struts 
centered  over  hulls  with  elliptical  cross  sections;  this  corresponded  to  the  configurations  of  interest  at  that  time. 
Subsequent  advances  in  computers  diminished  the  motivation  for  approximations  and  a  wider  variety  of  cross 
sections  was  considered  as  the  SWATH  concept  matured.  Consequently,  the  Frank  Close  Fit  Technique^  was 
incorporated  in  the  U.S.  Navy  code  and  was  optionally  utilized  to  calculate  the  added  mass  and  damping 
coefficients.  However,  since  utilization  of  the  velocity  potentials  for  the  exciting  forces  and  moments  required  a 
reorganization  of  the  program,  the  wave  exciting  forces  and  moments  continued  to  be  calculated  as  functions  of  the 
added  mass  and  damping  coefficients. 

The  eomputer  program  which  implemented  this  model  is  known  as  the  SWATH  Motions  Program 
(SWMP).  Predictions  from  SWMP  for  the  vertical  plane  agreed  well  with  model  scale  experimental  results  for  all 
headings  through  moderate  speeds,  but  predictions  for  the  transverse  plane  responses  were  not  reliable.  Fortunately, 
for  many  SWATH  configurations,  the  roll  natural  period  is  in  a  range  where  little  wave  energy  occurs,  resulting  in 
little  roll  response.  However,  the  need  for  a  more  accurate  predictions  for  the  transverse  plane  was  clear. 

SWMP94,  a  recent  modification  of  SWMP,  implements  the  mathematical  model  which  is  described  in  this 
report.  In  SWMP94,  exciting  forces  and  moments  are  a  function  of  velocity  potentials  which  are  calculated  using 
the  Frank  Close  Fit  Technique.  The  mathematical  model  for  the  effects  of  appendages  differs  somewhat  from  the 
one  in  SWMP.  Expressions  to  define  the  lift  and  drag  coefficients  for  the  transverse  plane  are  developed.  These 
coefficients  vary  with  the  geometry  of  the  ship  and  are  determined  within  the  program.  Corrections  are  made  to 
Lee's^  derivation  for  the  cross  flow  drag  and  lift  components  for  the  transverse  plane.  Also,  some  of  Lee's 
simplifications  are  removed.  Consequently,  the  vertical  velocities  in  the  cross  flow  drag  contributions  to  the  forces, 
as  well  as  new  damping  terms,  couple  the  transverse  and  vertical  planes  of  motion. 

MATHEMA'nCAL  MODEL 

In  modeling  the  motions  of  SWATH  ships  in  waves,  a  right-handed  coordinate  system  which  has  x  positive 
forward,  z  positive  up  is  used.  Responses  are  calculated  about  the  longitudinal  center  of  gravity  of  the  ship  at  the 
ship's  centerline  at  the  mean  calm  waterline.  The  standard  linear  equations  of  motion  for  a  rigid  body  with  harmonic 


*  Presented  by  Dalzell  in  Stevens  Institute  of  Technology  reports  of  limited  distribution. 
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exciting  forces  are: 

(M  -h  A22  %2  +  ^llkl  +  {^24  “  %A  +  ^24^4  +  ^25^6  +  ^26^6  =  ^2^ 

(M  +  A33  )^3  +  633^3  +  C33^3  +  A35^5  +  635^5  +  035^5  =  F3e 

(I4  +  A44  )^4  +  644^4  +  C44^4  +  ( A42  “  Mzcg  %2  +  ^42^2  +  ^46^6  +  ^46^6  =  ^4® 

(h  +  A55  )^5  +  B5545  +  C55^5  +  A53^3  +  653^3  +  €53^3  =  F5e 

(h  "*"  ^66  )^6  +  ^66^6  +  ^62^2  +  ^62^2  +  ^64^4  +  ^64^4  =  ^6^ 

where  M  is  the  mass  of  the  displaced  volume  of  the  ship;  I4, 15,  are  the  mass  moments  of  inertia  in  roll,  pitch,  and 

yaw;  zcg  is  the  z  coordinate  of  the  ship's  center  of  gravity;  and  cOg  is  the  wave  frequency  of  encounter  where 
2 

^  TT  Q 

(Oe  =C0 - UCOSP 

g 

(D  is  the  wave  frequency,  g  is  the  acceleration  due  to  gravity,  U  is  the  forward  speed  of  the  ship,  and  (1  is  the  heading 
of  the  ship  relative  to  the  wave,  with  P  =  180°  for  head  seas.  Ajk  Bjk,  and  Cjk  are  the  added  mass,  damping,  and 
restoring  coefficients  in  the  jth  mode  due  to  the  k^h  motion  where  j  or  k  =  2,  3,  4,  5,  and  6  for  sway,  heave,  roll, 
pitch,  and  yaw,  respectively.  Fj  and  are  the  wave  exciting  force  and  displacement  of  the  ship  in  the  j^h  mode 

where 

The  hydrostatic  restoring  forces  are  defined  by: 

C33  =  pg  J  tdx 

C35  =  C53=PgJxtdx 
C44  =  pgj  y^tdx  -  MgBG 
C55  =  pgJ  x^tdx  -  MgBG 

where  t  is  the  thickness  of  the  strut  and  BG  is  the  distance  between  the  center  of  buoyancy  and  the  center  of  gravity. 
Surge  is  approximated  by: 
m|i  = 

The  coefficients  Aj^,  Bjk,  and  Cjk  and  the  forces  Fj  have  contributions  due  to  potential  flow,  cross  flow 
drag,  and  lift  due  to  the  body  and  the  appendages.  That  is,  in  general, 

D  =  (D)pp  +  (E))bl+bcd  (^)a 

where  D  may  be  Aj^,  Bj^,  Cjk  or  Fj,:  and  PF,  BL,  BCD,  A  stand  for  potential  flow,  body  lift,  body  cross  flow  drag, 
appendage,  respectively.  The  formulations  for  the  various  components  will  be  given  in  this  report. 

Potential  Flow  Components 

The  potential  flow  added  mass  and  damping  coefficients  and  the  wave  exciting  forces  for  sway,  heave,  roll, 
pitch,  and  yaw  are  those  presented  by  Lee^  which  result  from  application  of  the  strip  theory  of  Salvesen,  Tuck,  and 
Faltinsen^  to  SWATH  ships.  Potential  flow  terms  are  given  in  Appendix  A. 
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Body  Lift  and  Cross  Flow  Drag  Components 

The  mathematical  model  for  the  body  lift  and  cross  flow  drag  terms  which  will  be  presented  here  is  an 
extension  of  the  one  developed  by  Lee^.  Some  errors  in  Lee's  derivation  of  the  cross  flow  drag  and  lift  terms  for  the 
transverse  plane  have  been  corrected.  Also,  the  derivation  presented  here  differs  from  Lee's  in  that  the  relative 
vertical  velocities  include  heave,  pitch,  and  roll.  Consequently,  there  is  a  coupling  between  the  vertical  and 
transverse  planes  due  to  the  relative  vertical  velocity  terms,  resulting  in  the  introduction  of  the  terms  B34,  B43,  B45, 
and  B54. 

The  basic  development  for  the  body  lift  and  cross  flow  drag  components  is  given  in  Appendix  B.  Body  lift 
contributions  for  the  vertical  plane  were  developed  in  a  semi-empirical  manner,  as  summarized  in  Appendix  C.  A 
more  direct  approach  is  used  for  the  transverse  plane,  as  given  in  Appendix  D.  An  approach  for  evaluating  the  cross 
flow  drag  coefficients  is  given  in  Appendix  E. 

Appendage  Cross  Flow  Drag  and  Lift  Components 

Development  of  fin  contributions  to  the  coefficients  and  forces  for  stabilizing  fins  is  given  in  reference  3. 
This  model  was  used  in  references  1,  5,  and  6.  In  reference  6,  an  expression  for  the  lift  curve  slope  of  the 
appendages,  including  a  modification  for  the  downwash  from  the  forward  fin  on  a  coplanar  aft  fin  was  introduced. 
In  unpublished  work,  Ernest  E.  Zarnick  and  Ralph  G.  Stahl  of  the  David  Taylor  Model  Basin  extended  Lee's  fin 
model  to  the  case  where  the  appendages  may  be  canted  with  respect  to  the  horizon.  In  the  model  presented  here, 
these  canted  effects  have  been  retained.  However,  the  appendage  mass  contributions  have  been  separated  from  the 
added  mass  terms;  forward  speed  terms  which  had  been  included  have  been  removed.  Expressions  for  these 
components  are  given  in  Appendix  F. 

Total  Components 

The  contributions  to  the  added  mass,  damping,  and  restoring  coefficients  and  wave  exciting  forces  as 
described  above  are  summed  and  presented  in  Appendix  G. 

COMPARISON  BETWEEN  EXPERIMENTAL  AND  NUMERICAL  PREDICnONS 

In  this  section,  responses  to  regular  waves  which  were  measured  during  model  scale  experiments  and  those 
which  were  calculated  using  the  mathematical  model  described  in  this  report  are  compared  in  order  to  check  the 
validity  of  the  mathematical  model.  Characteristics  of  the  hullforms  for  which  results  are  presented  are  given  in 
Table  1.  The  hullforms  are  listed  in  the  order  in  which  they  were  tested.  The  hullforms  have  a  variety  of 
characteristics.  These  hull  forms  were  investigated  over  a  long  period  of  time,  reflecting  a  shift  firom  exploration  to 
fully  realized  ship  designs. 
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Table  1 .  Some  characteristics  of  hullforms. 


6A 

6B 

6C 

Stretched 

SSP 

TAGOS- 
81 -design 

TAGOS- 
81 -heavy 

TAGOS- 

19 

TAGOS- 

23 

Displacement,  LT 

2900 

2900 

2900 

614 

3014 

3463 

3338 

5370 

LOA,  ft 

240.0 

240.0 

240.0 

149.1 

244.0 

244.0 

232.50 

279.00 

Draft,  ft 

26.67 

26.67 

26.67 

16.5 

22.11 

29.20 

24.75 

26.00 

SD,ft 

75.0 

75.0 

75.0 

40.0 

77.10 

77.0 

71.00 

72.00 

Awp.  ft 

2121.8 

2123.0 

1908.0 

725.49 

2151.1 

2151.1 

2685.5 

4354.0 

LCB,  ft 

116.3 

115.3 

113.9 

72.11 

122.08 

122.2 

107.90 

129.90 

LCF,  ft 

113.85 

111.06 

108.10 

80.02 

122.67 

122.7 

107.41 

130.64 

GmL.  ft 

20.0 

38.0 

45.0 

50.52 

26.1 

38.2 

27.80 

49.00 

Gmt.  ft 

11.06 

10.60 

10.94 

5.25 

8.30 

14.1 

9.20 

8.23 

KG,  ft 

34.0 

34.0 

34.0 

14.11 

28.8 

25.3 

31.50 

33.80 

KB,  ft 

10.6 

10.6 

10.6 

7.25 

8.37 

10.53 

10.30 

11.18 

ke.ft 

55.91 

58.79 

58.07 

40.7 

57.0 

58.7 

41.36 

65.90 

k<|),  ft 

39.60 

39.76 

38.91 

18.4 

40.7 

41.7 

29.0 

38.8 

Struts/Hull 

1 

1 

2 

2 

1 

1 

1 

1 

Length  at  WL,  ft 

172.3 

208.0 

197.5 

139.0 

199.48 

199.48 

190.00 

230.00 

Strut  Setback,  ft 

27.7 

7.1 

9.4 

7.87 

16.47 

16.47 

16.00 

19.00 

HorizA^ert  Diam. 

1.0 

1.0 

1.0 

1.0 

1.8 

1.8 

1.4 

1.25 

WL  to  Hull 

11.67 

11.67 

11.67 

5.58 

9 

16.09 

8.75 

7 

WL  to  midhull 

19.17 

19.17 

19.17 

11.15 

15.56 

22.64 

16.75 

16.5 

A  comparison  between  the  experimental  and  predicted  results  are  given  in  Figures  1  to  8.  The  experimental 
results  for  the  6A,  6B,  and  6C  configurations  are  presented  in  Reference  16.  The  experimental  results  for  the  other 
configurations  are  presented  in  David  Taylor  Model  Basin  reports  of  limited  distribution.*  Experiments  were  run  for 
a  range  of  speeds  and  headings  for  the  various  configurations.  Bow  and  quartering  waves  were  both  45  degrees 
from  beam  seas.  The  order  of  the  experimental  results  corresponds  to  the  order  in  which  the  experiments  occurred. 


*  The  experimental  results  for  the  Stretched  SSP  and  the  1981  T-AGOS  design  were  reported  by  A. 
Gersten.  The  results  for  the  T-AGOS  19  were  reported  by  R.  Thomas  Waters  and  David  D.  Hayden  and  the 
experimental  investigation  for  the  T-AGOS  23  was  reported  by  James  W.  Hickok,  Jr.  and  Ralph  Stahl. 
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There  are  differences  between  the  conditions  which  are  assumed  in  the  mathematical  model  and  those 
which  are  possible  during  experiments.  In  the  mathematical  development  it  is  assumed  that  speed  and  relative  wave 
heading  are  constant.  In  the  experiments  it  can  be  difficult  to  maintain  these  conditions,  particularly  in  oblique 
headings.  Models  are  self-propelled  and  are  tethered  to  the  carriage  which  moves  above  the  model.  Speed  and 
deflection  of  control  surfaces  are  controlled  manually.  While  application  of  advances  in  electronics  have  made  it 
easier  to  maintain  speed  and  heading,  the  dominant  factor  in  control  of  the  model  is  speed.  Since  lift  is  proportional 
to  speed  squared,  the  ability  to  control  a  model  (or  ship)  increases  significantly  with  speed. 

For  some  of  the  configurations,  measured  roll  is  presented  for  head  and  following  waves.  In  head  or 
following  waves,  if  the  waves  are  purely  two-dimensional  and  the  heading  is  precise,  a  symmetric  ship  will  have  no 
roll  response.  Consequently,  the  magnitude  of  roll  in  these  conditions  gives  a  measure  of  the  extent  to  which  the 
intended  relative  wave  heading  conditions  were  not  met. 

Another  difficulty  is  that  under  some  conditions  the  model  will  surge  or  significantly  deviate  from  its 
intended  course,  resulting  in  taught  tethers.  This  is  particularly  challenging  at  zero  speed.  One  innovation  which 
was  introduced  in  the  more  recent  tests  is  used  for  beam  seas  at  zero  speed.  In  this  case,  the  model  is  turned  90° 
relative  to  its  usual  orientation  to  the  carriage,  waves  are  generated  from  the  short  bank,  and  the  model  is  allowed  to 
drift  with  the  waves  while  the  carriage  moves  with  it.  This  alleviates  the  interference  of  the  tethers  with  the  model. 

SWATH  responses  have  been  shown  to  be  nonlinear^^  with  respect  to  wave  amplitude,  particularly  at  low 
speed.  The  mathematical  model  presented  in  this  report  models  this  behavior  with  the  introduction  of  the  cross  flow 
drag  components  and  application  of  equilinearization.  When  calculations  are  made,  a  value  for  wave  amplitude  is 
needed.  In  comparing  predictions  with  experimental  results  it  is  convenient  to  specify  wave  slope  (the  ratio  of  wave 
height  to  wave  length).  In  the  predictions  which  are  presented  here,  two  values  of  wave  slope  are  used  in  order  to 
bracket  the  values  used  in  the  model  experiments.  Consequently,  there  are  two  curves  of  predictions  on  each  plot. 

CONCLUSIONS 

A  frequency  domain  mathematical  model  for  the  responses  of  SWATH  ships  to  regular  waves  has  been 
presented.  The  model  requires  a  description  of  the  hullform  geometry  and  mass  characteristics.  Expressions  have 
been  developed  for  explicitly  evaluating  the  various  lift  and  cross  flow  drag  components.  These  utilize  experimental 
results  from  a  range  of  sources.  Correlation  with  model  experimental  results  is  excellent  for  a  variety  of  hull  forms, 
across  a  range  of  headings  and  speeds. 
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Figure  1  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  SWATH  6A. 
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Figure  1  (Continued) 
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Figure  1  (Continued) 
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Figure  2  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  SWATH  6B. 
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Figure  2  (Continued) 
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Figure  2  (Continued) 
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Figure  3  -  Comparison  of  Experimental  and  Predicted  Motion  Transfer  Functions 
for  two  different  wave  slopes  for  the  SWATH  6C. 
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Figure  3  (Continued) 
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Figure  3  (Continued) 
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Figure  4  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  Stretched  SSP. 
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Figure  5  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  1981  T-AGOS  at  the  design  draft. 
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Figure  6  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  1981  T-AGOS  at  the  heavy  draft. 
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Figure  7  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  T-AGOS  19. 
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Figure  8  -  Comparison  of  experimental  and  predicted  motion  transfer  functions 
for  two  different  wave  slopes  for  the  T-AGOS  23. 
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APPENDIX  A  -  POTENTIAL  FLOW  TERMS 


In  the  equations  which  follow,  a]±  and  are  the  two-dimensional  added  mass  and  damping  coefficients. 


In  the  following  three  equations,  integrals  are  taken  around  a  two-dimensional  cross  section: 

^22  "71“  ^22  [^2  ^2^1 
C(x) 


iCOp 


^23--— b33-p  J(|)3n3di 
^  ®  C(x) 

^2^-- - b24=p 


C(x) 


n2  and  n3  are  components  of  the  unit  vector  normal  to  the  body  surface  and  (t)^  is  the  complex  velocity  potential  for 
forced  oscillation  in  the  mode. 

Reference  3  gives  the  potential  flow  contributions  to  the  added  mass  and  damping.  Integrals  are  evaluated 
from  the  tail  to  the  nose  with  the  origin  at  the  longitudinal  center  of  gravity.  Contributions  from  both  hulls  are 
included  in  the  integrations. 


(A22)pF=Ja22(x)dx 

(B22  )pF  -  J  b22(x)dx 

(^24)??  —  J  ^24(^)dx 

(B24)pF=Jb24Wdx 

(A26)pF  =  f  xa22(x)dxH  ^B22 

J  ©e 

(B26)pf  =  Jxb22(x)dx-UA22 

(A33)pf  =  Ja33(x)dx 

(B33)pf  =  J*  b33(x)dx 

( A35  )PF  =  -  f  xa33  (x)dx  -  B33 

■>  ©e 

(B35)pf  =  -  J  xb33(x)dx  +  UA33 

(A42)pF=Ja24(x)dx 

(B42  )pF  =  J  b24  (x)dx 

(A44)pf=  Ja44(x)dx 

(B44)pF  =  J  b44(x)dx 

(A46)pf=  fxa24(x)dx  +  -^B24 

J  ®e 

(B46  )pF  =  J  xb24  (x)dx  -  UA24 

( A53  )PF  =  -  f  xa33  (x)dx  +  B33 

COg 

(B53)pf  =  -  J  xb33(x)dx  -  UA33 

f 

(A55)pf  =  J  x2a33(x)dx +— A33 

COg 

f  9 

(B55)pf  =  1  X  b33(x)dx  +  — yB33 

(A62)pf  =  f  xa22(x)dx-^B22 

•>  ©e 

(B62  )pf  =  J  xb22  (x)dx  +  UA22 

(A64)pf  =  J  xa24(x)dx--^B24 

(B64)pf  =  J  xb24(x)dx  +  UA24 

f  2 

(A66)pF  =  ^  a22dx+  2  ^22 

J  ©e 

f  2 

(®66)PF  =  1  ^  b22dx4  ^B22 
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With  symmetry  of  the  ship  about  the  ship's  centerline  assumed,  the  contributions  to  the  wave  exciting 
forces  due  to  potential  flow  are: 

(F2)pf  = -2pgA Jdxe*''’“'°®P  J  |in2sin(kysinP)  +  ^[n2sin(3cos(kysinP)  +  n3sin(kysinP)](t)2|e'‘''dl 


(F3)pp  =2pgAjdxe'’'’“^°®^  J  |n3Cos(kysinP)  +  i  — [n2sinPsin(kysinp)-n3COs(kysinP)](j)3|e'“dl 

C(x)  ® 

(F4)pf  =-2p8Aj<lxe"“'“l> 

f  |i[yn3-zn2]sin(kysin|J)+^[n2sinPcos(kysinP)+n3Siii(kysinP)]»t'4je*^Ml 

C(<)  “ 

(F5)pF=-2p8A|<lxe»“»l' 

I  |xn3Cos(kysmP)  +  i^  [nj siiiPsin(kysinP)-n3Cos(kysinP}](>3|e'^dl 


(F6)pF=-2pgA|dxe“““f 

f  Jixn2  sin(kysinP)  +  “  [n2sinPcos(kysinP)  +  n3sin(kysinP)](t)2  ie^^dl 

J  (0  iCOp  I 

rtvti  L  ej  } 


Approximations  for  the  surge  wave  exciting  force  and  surge  wave  exciting  force  due  to  pitch  are: 
(Fj)pp  =  -2ipgAkcospJ|J  cos(kysinP)dzdydx 

(F51  )pF  =  2ipgAk  cos  P  JJI  ^  cos(  ky  sin  P)dzdydx 
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APPENDIX  B  -  CROSS  FLOW  DRAG  AND  BODY  LIFT  TERMS 


Lift  and  Cross  Flow  Drag  Vertical  Plane  Force 

From  Thwaites^,  the  vertical  force,  Fz,  for  an  inclined  body  at  a,  a  moderate  angle  of  incidence  relative  to 
the  flow,  is: 

^  *2 

where  p  is  the  mass  density  of  the  fluid,  Ay  is  the  area  projected  in  the  vertical  plane,  U  is  the  forward  speed  of  the 
ship,  aov  is  the  lift  coefficient,  a  is  the  angle  of  incidence  of  flow,  C^y  is  the  cross  flow  drag  coefficient  and  w  is 
the  relative  vertical  fluid  velocity  with  respect  to  the  body.  From  Figure  B.l  it  follows  that 
sin{a-^5)  =  ^ 

W  =  -Zb 

so  that,  for  small  angles, 
w  =  U(a-^5) 


Z 


Figure  B.l .  A  sketch  of  a  body  moving  with  respect  to  the  fluid  in  the 
vertical  plane. 

For  an  arbitrary  point  (x,y,z), 

w  =  43-45+y44-Cv(x.y.z) 

where  the  vertical  velocity  of  the  fluid  induced  by  the  incoming  wave  is 
dz 

and  (|)o,  the  complex  velocity  potential  for  the  incident  wave,  is  defined  by 

X  _  igA  kz+ikxcosp-ikysinp 
<1,0— —e 

where  g  is  the  acceleration  due  to  gravity,  A  is  the  wave  amplitude,  O)  is  the  wave  frequency,  k  is  the  wave  number, 
P  is  the  heading  of  the  ship  relative  to  the  wave  and  (Oe  is  the  wave  encounter  frequency.  Consequently, 

w  =  43  -  x45  +  y44  +  p  (B  2) 

a  =  ^5  +  (43  -  x45  +  y44  +  )  /  U  (B.3) 

Equations  (B.2)  and  (B.3)  can  be  substituted  in  Equation  (B.l)  to  define  Fy. 
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T  jft  and  Cross  Flow  Drap  Transverse  Plane  Force 

Applying  Thwaites'"^  expression  to  the  transverse  plane,  the  horizontal  force  Fy  for  a  body  with  a  moderate 
drift  angle  is: 

Fy  =-^AH(u^aoHy  +  CDHv|v|j  (B-4) 

where  Ah  is  the  area  projected  in  the  horizontal  plane,  aoH  is  the  lift  coefficient,  yis  the  angle  of  incidence  of  flow, 
Cdh  is  the  cross  flow  drag  coefficient,  and  v  is  the  relative  fluid  velocity  with  respect  to  the  body.  Referring  to 
Figure  B.2, 

sin(-Y  +  ^6)  =  ^ 
v  =  -yb 

so  that,  for  small  angles, 

v  =  _u(y-^6) 


y 


Figure  B.2.  A  sketch  of  a  body  moving  with  respect  to  the  fluid  in  the 
transverse  plane. 


For  an  arbitrary  point  (x,y,z), 

V  =  42+46-44-CH(x.y.z) 

where  the  horizontal  velocity  of  the  fluid  induced  by  the  incoming  wave  is 


Ch  = 


90O  .-im.t 

dy 


Consequently, 

V  =  42  +  x46  -44  +coAslnpe''^+‘*'’‘=“P-'‘‘y*=“P 
7  =  ^6 -(^2+46-44  +toAsinPe''^^“'’‘^“P-'''y*‘"P)/U 

Equations  (B.5)  and  B.(6)  can  be  substituted  in  Equation  (B.4)  to  define  Fy. 


Lift  and  Cross  Flow  Drag  Moments 

The  roll ,  pitch,  and  yaw  moments  (Mx,  My,  Mz)  can  be  defined  as: 

Mx  =yFz+zFY 

My  =-xF2 
M2  =  xFy 


(B.5) 

(B.6) 
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(B.7) 

(B.8) 

(B.9) 


Evaluation  of  Body  Cross  Flow  Drag  and  Lift  Components 

Equations  (B.2)  and  (B.3)  are  substituted  in  equation  (B.l)  and  equations  (B.5)  and  (B.6)  are  substituted  in 
equation  (B.4)  and  applied  along  the  ship.  When  terms  are  collected  depending  on  whether  they  include  velocities, 
displacements,  or  terms  related  to  the  incident  wave  potential,  contributions  to  damping,  restoring  coefficients,  and 
wave  exciting  forces  are  identified. 

In  order  to  evaluate  the  forces  Fy,  F^,  Mx,  My  and  for  the  entire  ship,  the  terms  which  are  defined  in 
equations  (B.l)  to  (B.9)  are  integrated  along  the  port  and  starboard  hulls.  For  example. 

This  equation  is  expanded  with  the  lift  and  drag  coefficients  assumed  to  vary  longitudinally.  Terms 
associated  with  cross  flow  drag  coefficients  include  terms  of  the  form  b|b[,  where  b  is  a  velocity.  These  terms  are 

nonlinear;  consequently,  in  order  to  be  able  to  solve  these  terms  in  the  linear  domain,  equilinearization  is  applied. 
That  is,  for  any  harmonic  motion  given  by  b  =  bo  costOgt ,  b|b|  can  be  approximated  in  the  following  manner^: 

♦  1*1  8 

bb  = — (Ofibbo 
'  '  3n 

For  the  port  hull  Fz  is  applied  at  (x,  Sor  di)  and  for  the  starboard  hull  it  is  applied  at  (x,  -Sd»  -di)  where 
Sd  is  the  distance  from  the  ship  centerline  to  the  strut  centerline,  and  di  is  the  distance  between  the  calm  waterline 
and  the  maximum  beam  of  the  hull  or  half  the  draft  fpr  sections  with  no  lower  hull.  The  projected  area  Ay  in 
equation  (B.l)  is  replaced  by  integrating  dn,  the  horizontal  diameter  of  the  lower  hull.  When  equilinearization  is 
applied,  and  terms  are  grouped,  Fz  produces  the  following  contributions; 

(B33  )bL+BCD  ~  J  ^  ^  J  (CdVpIwp  |  +  Cpys  l^s  |)ix 

(B34  )bl+bcd  ~  f  1  (Cdvp  1  wp|  -  Cdvs  ks  |)ix 

(B35  )bl+bcd  =  -pU  J  xaovdndx  “I-— JxdH(CDVp|wp|+ Cdvs  ks  |)dx 

(C35)bL+BCD  =  pu^  JaovdH<^X 

(F3)bl+BCD  =  -ipUtoAj  aovdH  cos(kSD  sinp)e-''‘*'-"““=‘’*Pdx 

-icoA^-^  [ dH(cDVplwp|e--'‘^‘’*‘"P  +CDVs|ws|e“‘^"*“^)e"‘'‘‘'^‘'““''Pdx 
2  3ti  J  V  ' 

where  aov»  ^DVs  vertical  lift,  port  hull  vertical  cross  flow  drag  and  starboard  hull  cross  flow 

drag  coefficients.  The  relative  vertical  velocities  at  the  port  and  starboard  hulls  are: 

Wp  =  43  -x45  +Sd44  + 

Ws  =43  -x45  -Sd44  + 

Similarly, 

FY=-^J{FYp+FYddx 
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is  applied  at  (x,  Sq,  -62)  and  (x,  -Sd,  -d2)  for  the  port  and  starboard  hulls  where  d2  is  defined  to  be  equal  to  dj  for 
fully  submerged  sections  and  is  d/2  for  other  sections,  where  d  is  the  transverse  projection  of  the  body.  The 
projected  area  in  equation  (B.4)  is  replaced  by  integrating  d.  The  following  contributions  result: 

(^22 )bL+BCD  “PU Jd(CDHp KpI CdHsI''s|)^^ 

(B24)bL+BCD  =~puJ^2^0H^‘^’‘‘^'^'^J‘^2^(f'DHpKp|'''^DHsl''s|)^^ 

(^26  )bL+BCD  “  -pu  J  xagH^dx  J  ^^(Cdhp  l^p  I  +  CpHs  ks 

(C26  )bL+BCD  =  pu^  J" 

(F2)bl+bcd  =  pUmAsinpJ  aoHdcos(kSD  sinP)e“''‘'2'*'‘''*“®'^dx 

2  3iz  J  ^  ' 

where  apH  and  C^Hp  and  C^hs  are  the  horizontal  lift  and  horizontal  cross  flow  drag  coefficients,  and  the  relative 

velocities  for  the  port  and  starboard  hulls  are: 

vp  =^2  +  x46  +d244  +  coAsinpe-'‘‘^^+''^’‘"“P-‘''^“*“P 

Vs  =  ^2  +  46  +  dlL  +  ®Asin  ^^-^<i2+ikxcosfi+ikS,  sinp 
Utilizing  equation  (B.8)  for  My,  where  xc  is  defined  in  Appendix  D. 

(®53  )bl+bcd  =  “PU  J  xagy^H'^x  -  — J  xdjj  (c^vp  I  wp  I + Cdv  s  Ks  |)ix 

(B54  )bL+BCD  B54  “  “  2  ^  J"  P  kp  I  “  ^DVS  I'^'S  |)dx 

(®55)bL+BCD  -  puj X^aovtlH^X  +  ^  — Jx^dH(CDVp|'''p|  +  CDVsl'''s|)^X 

(C55)bL+BCD  ~p^^  J  xaov^fH^X 

(F5)bl+BCD  =  ipUcoAj  xaovdH  cos(kSD  sinP)e"‘'‘‘''^‘’“‘'°*Pdx 

J  xdH(cDVp|wp|e-‘'^"'>“"P  +CDVs|ws|e'""““"P)e-“--^‘^ 

Utilizing  equation  (B.9)  forMz, 

(®62  )bL+BCD  ^  ”P^  J  ^  ^  1  ^^(^DHp  hp  I  ^DHs  Ks 

(®64)bL+BCD  =  “P^J  J^^2^(^DHpKpK^DHsKs|)^^ 

(B66  )bl+BCD  =  -pu  J  xcxaoHddx  +  ^  ^  J  x^d(CDHp  |vp  |  +  Cdhs  I  Vs  |)dx 

(C66)bL+BCD  “  pU^  JxcaoH^^^X 

(F6)bl+BCD  =  pUcoAsinpJ  xcaoHdcos(kSD  sinP)e“‘^‘*^‘^'^‘^“Pdx 

P  8  .  ■  of  1  I  -ikSnSinp  1,  |  ikSnSinpU -kdj+ikxcosp . 

coAsinpJxd|CDHplvp|e  '^+Cdhp|vsF  f  dx 

For  lift  terms,  the  moment  arm  xc  has  been  used.  This  is  based  on  the  concept  that  since  the  transverse  projection  of 
the  SWATH  can  be  thought  of  as  a  wing,  the  moment  associated  with  lift  will  not  act  about  the  longitudinal  center 
of  gravity,  but  more  nearly  the  quarter  chord  of  the  strut.  The  definition  of  xc  is  given  in  Appendix  D. 
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From  equation  (B;7)  for  Mx,  with  the  moment  arm  for  Fy  being  -d2  and  the  one  for  Fz  being  ±Sd,  it 
follows  that: 

(B42  )bL+BCD  =  “Pu|  +2^1  ‘^2d(CDHp  Kp  I  +  CDHsl''s|)dx 

(B43  )bl+bcd  =  J  SodH  (Cdvp  I Wp  I  -  Cdvs  Ks  l)dx 

(B44)bl+BCD  =  pu|  S^aovdndx  +  puj  d^aonddx 

+  ^^|sDdH(CDVp|wp|  +  CDVs|ws|)dx  +  ^  — Jd2d(CDHp|vp|  +  CDHs|vs|)dx 
(B45)bL+BCD  =  ”  2^1  ^^DdH(CDVpKp|-CDVs|ws|)dx 
(B46  )bl+BCD  =  -pu  J  xd2aoHddx  +  ^  ^  J  xd2d(CDHp  |  vp  |  +  Cqhs  Ks  |)dx 
(C46  )bl+bcd  =  PU^  J"  d2aoHddx 

(F4)bl+BCD  ~  ~PU®AjSj)aovdH  sin(kSD  sinP)e  ^dx 

+pUG)Asin(iJd2aoHdcos(kSj) 

— ®AjSodH(^Covp|wp|e  ^  ^-CDVsrsF 

_  p  ±  a)A  sin  p  f  d2d(cDHp  I  vp  +  Cdhs  |  vs  )e-‘=‘*^  +ikxcosp^^ 

2  3n  J  ^  ' 


In  order  to  evaluate  these  terms,  values  for  the  lift  and  cross  flow  drag  coefficients  must  be  defined. 
Vertical  plane  body  lift  components,  horizontal  plane  body  lift  components,  and  cross  flow  drag  coefficients  are 
developed  in  Appendix  C,  D,  and  E,  respectively. 
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APPENDIX  C  -  VERTICAL  PLANE  BODY  LIFT  CONTRIBUTIONS 


In  Appendix  B,  the  vertical  plane  cross  flow  drag  and  body  lift  components  for  the  vertical  plane  include 
body  lift  and  cross  flow  drag  coefficients  aoH'  Experimental  results  indicated  that  a  consistent  definition  for  aoH  led 
to  invalid  results.  Consequently,  semi-empirical  expressions  developed  for  submarine  derivatives*  which 
correspond  to  the  damping  terms  were  utilized  as  the  basis  for  the  body  lift  terms.  The  adaptation  of  these 


expressions  to  SWATH  configurations  is  described  in  reference 


(®33)bL  --pL^u|2. 


6. 


(B35)bl  =(B53)bl  =0.207pL^U(k2  -ki)m’ 


(B55)bL  - 


(BbsIbl 

(B33)bL 


(C35)bL  -U(B33)bl 

(C55)bl  = -0-278pL3u^k2  -  ki)m' 

where  L  is  the  ship  length,  m’  is  the  displaced  volume  of  the  unappended  ship  divided  by  pL^  and 


For  the  wave  exciting  forces,  aov  should  be  defined  as  follows: 


aov  - 


(B33)bl 

pjdndx 


*  Developed  by  Elizabeth  Dempsey  of  the  David  Taylor  Model  Basin  and  reported  in  a  report  of  higher 
classification. 
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APPENDIX  D  -  TRANSVERSE  PLANE  BODY  LIFT  COEFFICIENTS 


Transverse  Plane  Body  Lift  Coefficients 

The  portion  of  a  SWATH  configuration  which  is  below  the  waterline  can  be  described  as  being  composed 
of  sections  which  have  both  a  lower  hull  and  strut,  only  a  submerged  lower  hull  or  no  lower  hull.  In  the 
development  below,  sections  which  do  not  have  a  strut  are  assumed  to  be  either  nose,  tail,  or  parallel  sections.  In 
defining  lift  coefficients,  a  component  approach  has  been  taken.  Each  sectional  component  is  assumed  to  have  a 
constant  value  for  the  lift  coefficients  and  is  taken  to  be  generating  a  moment  about  the  section's  center  of  pressure. 

Nose  Section 

Utilizing  Figure  11  in  Hoerner^,  page  19-8,  the  lift  coefficient,  aoH»  for  a  plain  round  fuselage  can  be 
estimated  as: 

aoH  =  0-009Lsection  /  ^section 

where  LsecTION  ^SECTION  the  length  and  maximum  diameter  of  the  section.  The  data  is  for  sections  with 
LseCTION®SECTION  ranging  from  about  6  to  9.  From  Reference  8,  page  19-17,  the  center  of  lift,  Xc,  relative  to 
the  nose  of  the  section  is: 

Xc  =0.4Lsection 


Parallel  Section 

The  lift  will  be  neglected  since  it  is  theoretically  0.0,  and  will  be  small  compared  to  other  contributions. 

Tail  Section 

Based  on  Reference  8,  page  19-2: 

Uqjj  =  0.022 

From  Pitts  et  al^.  Chart  9,  the  center  of  pressure  of  an  ogival  nose  from  the  maximum  diameter  of  the  tail  is 
approximately, 

Xc  =0.55Lsection 

Strut  Section 

If  the  strut  section  is  taken  to  be  a  lifting  surface,  then,  from  Whicker-Fehlner^^  for  a  low  aspect  ratio 

control  surface  in  the  free  stream, 

1.  SrcARg 

aoH  = - /  7  . . 

1. 8  +  'Y 

For  the  SWATH  configuration,  ARg  is  the  aspect  ratio  of  the  strut,  defined  as  the  average  draft  divided  by  the  length 
of  the  section.  The  center  of  pressure  is 
Xc  =0-25Lsection 
from  the  leading  edge  of  the  section. 
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APPENDIX  E  -  CROSS  FLOW  DRAG  COEFFICIENTS 


Covp.  Cdvs>  Cdhp’  appear  in  the  expressions  for  the  cross  flow  drag  contributions  to  the 

coefficients  and  forces  which  are  presented  in  Appendix  B.  Two  sets  of  data  are  used  to  evaluate  these  coefficients. 
The  first  set  of  data  is  for  oscillating  plates  presented  by  Bearman  and  Graham' '  and  is  given  as  a  function  of  the 
Keulegan-Carpenter  number,  KC.  A  tabulation  of  these  data  is  given  in  Table  E.l.  The  second  set  is  for  oscillating 
cylinders  reported  by  Sarpkaya'2  is  tabulated  in  Table  E.2.  Sarpkaya's  data  is  presented  as  a  function  of  KC, 

and  a  frequency  parameter,  PkC>  where 

V  T 
KC^-li!!32LL 

d 


where  V^ax  is  the  amplitude  of  the  velocity,  T  is  the  period  of  oscillation,  D  is  the  diameter  and  v  is  the  kinematic 
viscosity  of  the  fluid. 

The  cross  flow  drag  coefficient  is  assumed  to  vary  along  the  length  of  the  ship.  In  applying  the  data,  Vmax 
is  the  amplitude  of  the  relative  velocity  ,  Vp,  Vg,  Wp,  or  Wj,  as  appropriate.  T  is  defined  using  the  wave  frequency  of 
encounter.  For  evaluation  of  C^vp  and  C^vs.  the  data  for  an  oscillating  cylinder  is  used  with  D  being  the 

horizontal  diameter  of  the  lower  hull,  dn-  These  data  are  modified  by  an  expression  to  reflect  the  influence  of  the 


strut  on  the  coefficient^. 


Cdv  - 


1.0- 


td^ 


Cd(U^,T,D,v) 


where  t  is  the  horizontal  thickness  of  the  strut,  dy  is  the  vertical  diameter  of  the  lower  hull. 

Cdhp  andCoHsare  evaluated  using  either  the  plate  or  the  cylinder  data.  For  sections  where  a  strut  is 

present,  the  plate  data  is  used  with  the  local  draft  used  for  D.  For  sections  without  a  strut,  the  oscillating  cylinder  is 
used  with  D  defined  as  dy. 


Table  E.  1 .  Cross  flow  drag  coefficient,  Cd  as  a  function  of  the  Keulegan  Carpenter  number,  KC,  for  flat  plates, 

based  on  data  presented  in  reference  11. 


KC 

Cd 

2.0 

6.75 

3.0 

5.75 

4.0 

5.15 

5.0 

4.75 

7.5 

3.90 

10.0 

3.40 

15.0 

2.95 

20.0 

2.80 
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Table  E.2.  Drag  coefficients  as  a  Function  of  Keulegan  Carpenter  number  and  a  frequency  parameter,  based  on  data 

presented  in  reference  12. 


Pkc  I 

KC 

497 

1107 

1985 

3123 

5260 

3.0 

1.64 

1.40 

na 

0.85 

5.0 

1.80 

1.65 

Is 

0.95 

10.0 

2.20 

1.95 

1.70 

1.45 

1.00 

12.5 

2.40 

2.05 

1.90 

1.30 

0.90 

15.0 

2.20 

2.00 

1.65 

1.10 

0.85 

20.0 

2.00 

1.70 

1.20 

0.90 

0.75 

30.0 

1.70 

1.40 

1.00 

0.70 

0.60 

60.0 

1.50 

1.05 

0.70 

0.55 

0.50 

150.0 

1.20 

0.90 

0.50 

0.48 

0.42 
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APPENDIX  F  -  APPENDAGE  CONTRIBUTIONS 


The  contributions  of  the  appendage  to  the  added  mass,  damping,  and  restoring  coefficients  and  wave 
exciting  forces  are  given  in  this  appendix.  Summations  are  taken  over  all  appendages. 

Let  cj  and  bi  be  the  mean  chord  and  span  of  the  appendage;  let  Ai  be  the  sweep  angle  of  the  quarter  chord 
line;  let  Xf  ,  yf  ,  Zf  be  the  coordinates  of  the  point  on  the  appendage  which  is  the  intersection  of  the  quarter  chord 

and  the  mean  span;  let  aj  be  the  cant  angle  of  the  appendage  relative  to  the  horizon.  Assuming  an  elliptical  cross 

section,  the  area  and  added  mass  can  be  defined  as: 

Af.  =Cibi 

A 

af  =— CjAf 

is  evaluated  using  the  data  given  in  Appendix  E  for  oscillating  plates.  While  alternate  values  for  the 

lift  curve  slope,  Clo.  .  can  be  specified,  the  default  expression  for  stabilizers  and  rudder  is: 

^  _  I.SjcAEj _ ^ 

“  j  ^ 

AR- 

1.8  +  cosAi- - —  +  4.0 

Vcos^Ai 

If  Clc  =  1.  this  is  the  expression  developed  by  Whicker  and  Fehlner^O  and  utilized  by  Dempsey  for 
submarines.  jg  the  effective  aspect  ratio  of  the  appendage,  which  is  based  on  the  tip  to  tip  span  of  a  pair  of 

appendages.  For  SWATHs,  unlike  submarines,  appendages  are  not  used  in  pairs.  However,  comparison  with 
experimental  results  showed  that  evaluation  of  CLcii  assuming  a  pair  of  fins  while  evaluating  the  lift  based  on  the 

actual  stabilizer  worked  well  for  SWATHs.^  Clc  is  a  modification  factor^  to  reflect  the  effect  of  downwash  from 
one  appendage  on  another  in  the  same  plane  -  that  is,  of  one  stabilizer  on  another.  Therefore,  Clc  equals  1.0  except 
when  the  appendage  is  aft  of  a  coplanar  appendage.  Data  presented  by  Lloyd^^  was  corrected  for  boundary  layer 
effects  and  tabulated  by  Cox  and  Lloyd.^^  Clc  is  given  in  Table  F.l. 


Table  F.l  Ratio  of  lift  on  aft  fin  to  lift  on  forward  fin  for  various  fin  separations  and  oscillation  frequency-to-speed 

ratios  (from  references  14  and  15). 


x/b 

rob/U 

0.00 

0.04 

0.08 

0.12 

0.16 

0.20 

10 

0.412 

1.221 

15 

0.462 

0.638 

0.846 

1.046 

20 

0.529 

0.732 

1.000 

1.151 

1.110 

0.971 

25 

0.614 

0.816 

1.099 

1.132 

1.101 

0.897 

30 

0.706 

0.853 

1.118 

1.006 

0.912 

0.853 
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(A22)A=X®f. 

(•^24 )a  =  -^,af  sinaj^Zf.  sin  a;  +yf.  cosajj 

(A26)A  =  XH^fi  “i 

(A33)a  =X^fi  “i 

(A35)a 

(A42)a  =-y^.af  sinai|zf.  sinaj  +yf,  cosaj  j 

(A44)a  =^af.  (zf,  sinaj+yf.  cosa;)^ 

(A46)a  =-^Xf  af  sinaj^Zf.  sina;  +yf.  cosajj 

(A53)a 

(A55)A=X^fi^fi^°®^  “i 
(A62)A=X^f|^fi 

(A64)a  =-y^,Xf  af_  sinttj^Zf.  sin  a,  +yf,  cosajj 
(A66)a  “i 

(B22)a  =f  sin^  aj 

(B24)a  =-|u]^Af,CLa,  sinai(zf.  sintt;  +yf.  cosaj) 

(B26)a  =|uJ^Xf,Af.CLa.  sin^aj 

(B33)a  =f  U^Af  Cl_^^  cos^tti  +£i- J^Af^  cosai(cD^  K|+Cd,^  |wpj) 
(B35)a  =-f  ‘=0®“i(^DjwsJ  +  CD_Jwp^ 

(B42)a  =-|u]^Af.CL„.  sinai(zf,  sina;  +yf^  cosaj) 

(B44)a  =|u^Af  CL„.(zf.  sintti  +yf^  cosaj)^ 

(B46)a  =-|u5],Xf.Af,CLa,  sinai(zf.  sinaj  +  yf.  costti) 

(B53)a  cosai(cDjw,J+CD^^  |wp^ 

(B55)a  =f  cos^a;  +£^]^x?  Af^  cosai(cD^^  KpJ) 

(B62)a  “ Af.CLg^.  sin  (Xj 

(B64)a  =-|u2],Xf. Af  CLa.  sinai{zf.  sintt;  +yf.  cosaj) 

(B66)a  =|u^x?  Af  Clk.  sin^tti 


(C26)a  sin^a; 

(C35)A=-|u2X^f.^L„, 

(C55)a  =-^u22^Xf  Af  Cl^^  cos^tti 

(C46)a  =-^u2^Af  CLot^  sinai(zf^  sinaj  +yf^  cosaj) 

(C66)a  =f  sin^aj 

(F2)a  =-AU(O^Af,CL(x.  sinaj^cosaj sin(kyf,  sin  p)  -  sin  p  sin  a;  cos{kyf  ^  sinpjp^'* 

costti  cos(kyf,  sinp)  +  sinpsinaj  sin^kyf.  sinpjje'“'' 


(F3)A=-ifAUco2;Af  C 


COS  a 


-i  — —  Acoy  Af  COStt: 
2  3n  ^  ‘ 


Wsje 


ik(xf.  cosp+yf.  sinp) 


+  Cd^^  wpje 


ik(xf.  cosp-yf.  sinpj 


-kZf 

k  M 


(F4)a  =-|AUcoJ]Af  CLa.  (zf,  sina;  +yf,  cosaj) 

cosajsin^kyf,  sinpj-sinPsinaj cos^kyf,  sinpjje*'^^'^**™'* 

(F5)a  =i-AU(O^Xf.Af,CL„,  cosai|cosai  cos(kyf.  sin  p)  +  sin  p  sin  aj  sin(ky f _  sin  p)]e‘'"f‘  ' 


+i— — Awy  Xf  Af 
2  371  ^  ' 


cosaj 


ik(xf.  cosP+y.,  sinp]  I  , 

J  V  f.  "^^  +  €0^  jWpJe 


ik(x,|  cosP-y(|  sinp)'\  -kz,. 


where 


(F6)a  =  2^Uco^Xf.Af  Cl„,  sin«i  cosaj  sin(kyf,  sinpj-sinPsinaj  cos(kyf,  sinp)je'“'' 

■  •  •  .  -kZf +ik(xf  cosp-Vf  sinp) 

wp,  =^3-Xf,^5+yf,^4+i®e  ■  ^  ^ 

•  •  t  •  -kZf +ik(xf  cosP+yf  sinp) 

Ws,  =^3-Xf,^5-yf,^4+i®e  ■ 


+ikXf.  cosp 


In  the  derivation  it  has  been  assumed  that  the  appendages  have  been  included  in  the  mass  and  inertias  of  the 

ship  and  that  the  origin  of  the  coordinate  system  is  at  the  appended  ship's  longitudinal  center  of  gravity.  In  general, 

the  effect  of  the  fin  on  the  mass  matrix  will  be  small.  However,  in  the  early  design  phase  the  mass  properties  of  a 

ship  is  estimated  and  under  some  circumstances  the  appendages  may  be  relatively  important.  In  that  case,  the  mass 

of  the  appendages  can  be  estimated  by: 
pit  . 

The  contributions  of  the  fins  to  the  mass  matrix  can  be  estimated  by: 

(1^22  )a 
(m24)A 

(m26)A 

(1^33 )a  - 
(11135)  A 
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(m44)A=2H'^yf.h^ 

(0146  )a 

(m53)A=-Xxfinif, 

(m55)A=XK 

(m62)A=2H"'fi 

(m64)A=-X^n^fi‘"fi 

(m66)A=X'"fi 


APPENDIX  G  -  TOTAL  COEEFICIENTS  AND  FORCES 


As  noted  in  the  main  text  of  this  report,  the  added  mass,  damping,  and  restoring  coefficients  and  wave 
exciting  forces  and  moments  are  modeled  with  potential  flow,  body  lift,  body  cross  flow  drag,  appendage  lift  and 

appendage  cross  flow  drag  components.  Details  of  the  development  of  these  components  are  given  in  Appendices 
A,  B,  C,  D,  E  and  F.  Expressions  for  the  evaluation  of  the  various  lift  (aoH.  aov-  ^Laj )  and  cross  flow  drag  (Cdh- 

Cdv.  )  coefficients  are  given  in  the  appendices.  Relative  vertical  velocities  are  defined  by: 

Vp  =^2  +46  +d244  +C0AsinPe-*^‘'^+‘'‘’‘=°^P-‘^"“*™P 
vs  =  ^2  +  46  +  ^2^4  +  COAsin  Pe-*^**^  +ikxcosP+ikS,  sinp 

Wp  =  ^3  -x^s  +Sd44  + 
ws  =  ^3  -x^5  -Sd44 

wp^  =^3-Xf,^5  +  yf.^4  +  itoe  ^  ‘  ’ 

•  •  •  -kZf +ik(x,  cosP+Vf  sinp) 

Ws,  =^3-Xf,^5-yf,^4  +  i®e  '  '■ 


The  totals  of  the  components  are  listed  below. 

A22  =  Ja22(x)dx  +  ^af.  sin^  aj 

A24  =  Ja24(x)dx-^af,  sinaj^Zf,  sinaj  +yf,  cosajj 
A26  =  Jxa22(x)dx+^Jb22(x)dx  +  ^Xf.af,  sin^Uj 
A33  =  J a33(x)dx  +  ^af,  cos^  aj 
A35  =  -  j  xa33  (x)dx  -  J  b33  (x)dx  "  X  “i 

A42  =  Ja24(x)dx-^af.  sinaj^Zf.  sinaj  +yf.  cosajj 
A44  =  I  a44  (x)dx  +  ^  af .  (zf ,  sin  ttj  +  yf .  cos  a; ) 

A46  =  Jxa24(x)dx  +  -^Jb24(x)dx-X,Xfiaf,  sinai(zf.  sinaj  +yf,  cosaj) 

A53  =  -  J xa33(x)dx  +  -^  Jb33(x)dx  -^Xf.af^  cos^  a; 

c 

2 

A55  =  J  x2a33  (x)dx  +  •;^  Ja33  (x)dx  +  X  “i 

A62  =  Jxa22(x)dx-^Jb22(x)dx  +  ^Xf^af^  sin^a; 

A64  =  Jxa24(x)dx--^jb24(x)dx-]^xf  af.  sinai(zf.  sinaj+yf^  cosaj) 

2 

Agg  =  J  x2a22dx  +  — I  a22  (x)dx  +  X?.  af .  sin^  aj 
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B22  =  Jb22(x)dx  +  ^U^Af,CLa,  sin^ 

-pU  J  aonddx  + I  d(Ci3Hp  I  vp  I  +  Cdhs  |''s  |)dx 

B24  =  Jl>24Wdx-^U^Af,CLa,  sinaj^Zf.  sinaj  +yf.  cosaj  j 
-puJd2aoHddx  +  ^^Jd2d(CDHp|vp|  +  CDHsl''s|)dx 
®26  =  Jxb22(x)dx-uJa22(x)dx  +  ^U^Xf  Af  Clo(.  sin^  a; 

-pu|xaoHddx  +  -^-^|xd(CDHp|vp|  +  CDHs|vs|)dx 

B33  =  J b33(x)dx  +  -^U^ Af  cos^  tti  cosaijCo^,  |+Cd^.  |wp^  |) 

-pL^ U^2. 439'’^^(k2  -kj)m’  - m' j  +  ^^jdn (c^vplwp |  +  CdvsIws |)dx 

B34  =  2  ^|^D‘^H(CDVp|wp|-CDVsl'''s|)dx 

B35  =  xb33  (x)dx  +  U  j  a33  (x)dx 

- f  Xf .  Af  cos^  ai  -  H ± ^  Xf ^  Af ^  cosai (Cd^^  |w,_  |  +  |wp^  |) 

+0.207pL^U(k2  -ki)m’--^— JxdH(CDVp|wp|+CDVsl'''s|)dx 
B42  =  j  b24  (x)dx  -  ^  Af  CLa,  sin  ttj  (zf .  sin  ttj  +  y  f ^  cos  tt; ) 

-pU  J  d2aoHddx  +  2  ^  J  *^2d(CDHp|vp  |  +  Cdhs^s  l)dx 

B44  =  Jb44(x)dx  +  |u2^Af.CL„,(zf.  sintti  +yf.  costti)^ 

+pU  J  SpaovdHdx  +  pU  J  d^aonddx 

+|^JsDdH(CDVp|wp|  +  CDVsK|)dx  +  |^jd^d(CDHp|vp|  +  CDHs|vs|)dx 

B45  =  J  xSDdH(CDVp|wp|-CDVs|ws|)dx 

B46  =  I xb24(x)dx-uj a24(x)dx-|u^Xf.  Af  CLa.  sinaj^Zf  sinaj  +yf.  cosajj 

-puj  xd2aoHddx  +  ^^  J  xd2d(CDHp|vp|  +  CDHshs|)dx 
B53  =  -  j  xb33(x)dx  -  uj  a33(x)dx 

-f  Af  “®“i(<^DjwsJ  +  CDjwpJ) 

40.207pL3u(k2-ki)m'-||^JxdH(CDVp|wp|  +  CDVs|wsl)dx 

B54  =  -^^JsDXdH(CDVplwp|-CDVsks|)dx 
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B55=Jx^b33(x)dx  +  ^jb33(x)dx 

X?  Af^  cos^  a;  +  £  a  ^  X?  Af  ^  cos  ttj  (Cd^^  |  +  |w 

-pL'^U — ^  J  ^  *1  ^ - +  ^^jx^dH(CDVp|wp|  +  CDVs|ws|)dx 

2.439'^^-^(k2 -kj)m' -m  2  3nJ 

B62  =  Jxb22(x)dx  +  uja22(x)dx  +  ^U^Xf,Af,CLa,  sin^ ttj 

-puJxcaoHddx  +  -^^Jxd(CDHp|vp|  +  CDHsl''s|)dx 


B64  =  J  xb24 (x)dx  +  u J  a24 (x)dx  -  Xf ^  Af ,  C^x,  sin tt;  (zf^  sin  ttj  +  yf ,  cosaj ) 

-puj  xcd2aoHddx+^— J'xd2d(Cp)Hp|vp|  +  CDHsl''s|)^x 
2 

-pujxcxaoHddx  +  ^^Jx2d(CDHp|vp|  +  CDHs|vsi)dx 


^26  + 

C33  =pgjtdx 

€35=  pg J xtdx  +  j X  i cos^  ttj  -  pL^U^ ^2. 439^^-^(k2-ki)m'  -  m'  j 
C44  =  pgj  y^tdx  -  MgBG 

C46  =  -^U^X^fi^Lai  sinai(zf,  sinaj  +yf.  cosajj  +  pU^  Jd2aoHddx 
C53=PgJ’‘tdx 

C55  =  pg  J  x^tdx  -  MgBG "J  ^  Xf .  Af .  cos^  ttj  •”  0. 278pL^U^  (k2  -  kj  )m' 

C66  =  y  ^  X?,  Af ^  Cl__,  sin^  ttj  +  pU^  J  xcaonddx 


With  symmetry  of  the  ship  about  the  ship’s  centerline  assumed,  the  contributions  to  the  wave  exciting 


forces  due  to  potential  flow  are: 

F2=-2pgAjdxe““‘=“P  J  I  in2  sin(kysin  p)  +  ^[n2  sinPcos(kysin  p)+  n3sin(kysinp)]<|)2|e'“dl 


C(x) 


+^AUto^Af  Clu.  sina;  cosa;  sin(kyf.  sinpj-sinPsinai  cos^kyf,  sinpjje*'^'*^^'* 
+pUa)Asin  P  J  aged  cos(kSD  sin  p)e“''‘*2''’^’‘®°®^dx 

_£  AcoAsinp  f  de-*'‘'^+“'"'“P(cDHp|vp|e-“‘®‘’*‘"^  +CDHsl'^s|e“'^"'“^)dx 
2  371:  J  '  ^ 
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F3  =  2pgAjdxe‘’^^“®^  J  |n3  cos(kysinP)  +  i  — [n2 sinpsin(kysinp)- n3 cos(kysinP)](j)3 je'^^dl 
C(x)  ® 

-i— AUco^Af  Cl^,  costti  cosaj  cos|kyf,  sinpj  +  sinpsinaj  sin^kyf,  sin 


2 

:P  8 


P^AcoYAf 

2  3ji  ^ 


( 


cosaj 


Cd,,  Ws,  e 

Sj  I 


ik(xf.  cosP+yf,  sinp) 


+  wpje 


ik(xf.  cosp-yf.  sinpj 


cosp 


-kZf. 


-ipUcoAj aov^H cos(kSD sinP)e  W,+ikxcosPjj^ 

_icoA^AjdH(cDVp|wp|e-‘''^'^“"P  +  CDVs|ws|e“'®'’™P)e-'‘^'-^‘'““*Pdx 


2  3ji 

F4=-2pgAjdxe'*'='^°*P 


J  I  i[yn3  -  zn2  ]sin(ky  sin  P)  + — [n2  sin  P cos(ky  sin  P)  +  n3  sin(ky  sin  P)](t)4  le'^^dl 

C(x)  ® 

-^AUco^Af  Clc  |zf.  sinaj  +  yf.  cosajj 

^cosa;  sin^kyf.  sinpj-sinPsinaj  cos|kyf,  sin pjje*'^'' 

-PUcoaJ  SpaovdH  sin(kSD  sinP)e”''‘*''‘’‘’“‘^°®^dx 

+pUcoAsinpJ  d2aoHdcos(kSD  sin  p)e“'^‘*^  “*^dx 

-i|^(oAjsDdH(cDVp|wple-‘‘'^'>“"P-CDvs|ws|e‘""°""P>-'‘‘''^‘'“ 


— ^^(0AsinpJd2d|cDHp|vp|e  +  CDHs|vs|e'*'^°*‘"^)® 

F5=-2pgAjdxe’'“"°*P 


-kdo+ikxcosf 


J  I  xn3  cos(ky  sin  [})  +  i  — 


U 

icOp 


[n2  sinPsin(kysinp)“-n3  cos(kysin 


C(x)  ' 

+i^AUa)^Xf.  Af.CLa.  cosai|cosai  cos^kyf.  sinpj  +  sinPsinaj  sin^kyf. 
•  P  8 


+i— — Ao)^  Xf  Af  cos  a; 
2  3n  ^  ' 


ik(x,.  cosP+y,.  sinp) 


+  Cd  Wp  e 

Pi 


ik(x,.  cosP-y,.  sinp) 


kz, 

*1 


+ipUo)Aj  xaoydn  cos(kSD  sin  p)e~''‘*'^^*“*^dx 

+i(0A^-^  J  xdHlCovpIwple-^^S-^^P  +CDVs|ws|e‘‘‘^‘’'*'"^)e-''‘*'-^"'’“=“Pdx 
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=  -2pgA  j  c 


^ikx cosP 


[  iixn2sin(kysinp)  +  —  x-^ —  [n2sinpcos(kysinP)  +  n3sin(kysinP)](t)2 
J  CO  iCOp 

C(x)^  L 

+^AUco^Xf,Af,CLaj  sinaj^cosaj  sin(kyf.  sinpj-sinPsinaj  cos(kyf^  ' 

+pUcoAsinpJ  xcaoH'^cos(kSD  sin  P)e“'‘‘’2  “®^dx 


_£i_a)AsinP  f  xd(cDHp|vp|e-'"'‘>""P  +CDHp|vs|e“^"'’""P)e-’'‘‘^-^“'"“*Pdx 
2  37C  J  '  ' 

Fj  =  -2ipgAkcosp|JJ  cos(kysinp)d2dydx 

F51  =  2ipgAk  cos  P  JJJ  cos(ky  sin  P)dzdydx 
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